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ABSTRACT 


*  The  equations  of  motion  of  nn  alrthlp  encountering  discrete 
una  rchuuiti  ya»:»  oirj  «.«►*<•  mu*.—  f  *  *L~  >■. u<-mrai  loads  are  fcumu- 

iated.  Tn6  aerodynamic  forces  and  moments  due  to  the  gust  and  the  resulting 
motion  are  derlvea  mom  »;«nder-uo dy  theory.  Corrections  are  applied  to 
account  for  viscous  effects.  A  numerical  example  Is  pre*ented  to  determine 
the  response  and  loads  of  a  typical  alrshio  penetrating  a  discrete  one-mtnus- 
coslne  gust  and  random:  turbulence.  The  discrete  gust  results  Indicate  that 
critical  loads  on  the  en.*alope  and  tall  are  cauied  by  a  gust  length  equoi  to 
the  length  of  the  airship.  The  random  gust  results  Indicate  that  the  critical 
scale  of  turbulence  Is  also  approximately  one  airship  length-  On  tho  hor.is 
of  these  computations  a  rational  procedure  In  airship  gust  load  cri'erio  is 
proposed. 
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CHAPTER  I 
INTRODUCTION 

1,1  Historical  Background 

Previous  work  on  gust  load  analysis  of  airships  is  liavtcd  and 
scattered.  The  earliest  significant  investigarion  in  this  field  was  conducted 
by  Munk  (Refs.  1  and!!)  who  mod  potent  lal'theory  to  calculate  rhe  forces 
and  moments  on  airship  hu1  Is  in  stoody  flow,  Further  theoretical  treatments 
using  , son-viscous  theories  were  pursued  by  Zahm  (Ref.  3)  ond  Upson  and 
Kli<off  (Ref.  4)  who  obtained  essentially  tho  some  results  os  Munk.  Some 
early  applications  of  Munk1*  potential  theory  in  determining  the  loads  im¬ 
posed  by  gusts  on  rigid  airships  ore  presented  in  Refs.  5  through  7. 

An  early  attack  on  the  origin  of  ilic  fores*  acting  on  dy  cf  ( 
revolution  was  made  by  Harrington  in  Ref,  jl.  Harrington  determined  experi¬ 
mentally  the  distribution  cf  vortlclty  in  the  wake  of  the  body  from  which  he 
wos  able  to  calculate  the  lift  by  using  Prandtl’s  lifting-line  theory.  His 
results  compared  fovorobly  wirh  measured  values.  In  1938  Kuethe(Pef.  9) 
used  a  wafer  tank  to  determine  the  forces  ond  moments  on  a  model  of  the 
Akron-Maoon  class  from  the  recorded  lateral  and  pitching  accelerations  of 
the  center  of  buoyancy.  A  gust  with  a  velocity  gradient  wos  simulated  by  a 
flow  in  a  channel  normal  to  the  path  of  the  model. 

Recent  investigations  into  thenaturo  of  the  forces  acting  on  bodies 
of  revolution  were  conducted  by  Alien  and  Perkins  (Ref.  10)  and  refined  by 
Kelly  (Ref.  11)  who  determined  *he  normal  force  and  pitching  moment  for 
angles  of  attack  beyond  the  range  of  potential  theory.  Their  method  con¬ 
sidered  the  additional  force  generated  by  the  viscosliy  of  the  cross-flow  and 
their  results  agreed  well  with  experiment.  At  about  the  same  time  Hill  (Ref.  17) 


ASRl  TR  72-1 


developed  a  theory  by  which  the  loading  was  predicted  by  replacing  the 
body  surface  with  o  vortex  sheet, 

Investigations  of  the  lift  and  moment  acting  on  a  slender  body 
performing  unsteady  motion  wure  conducted  both  theoretically  and  experi- 
mrntally  by  Ashley,  Zartarian  and  Neilson  (Ref.  13).  Foss  (Ref.  14)  deter¬ 
mined  expressions  for  the  growth  of  lift  and  moment  acting  on  a  slender  body 
as  it  penetrot»s  a  discrete  gust. 

mve-t 'gation  of  airship  gust  loadino 

was  conducted  by  Flomenhoft  (Ref.  15)  who  determined  the  transient  motion 
of  an  airship  in  response  to  discrete  gusts  striking  the  fins. 

1.2  Object  jf  the  Presenr  Investigation 

Current  design  methods  for  present  day  airships  are  based 
primarily  on  semi-empirical  techniques.  The  envelope  and  fins  are  designed 
using  the  concept  of  a  sharp-edged  gust  applied  directly  to  the  stern  of  the 
airship  (Ref.  16).  A  more  realistic  analysis  considering  the  transient  effect* 
resulting  from  the  groduat  penetration  of  the  envelope  into  a  gust  front  5s 
not  utilized. 

!t  is  the  purpose  of  the  present  investigation  to  apply  rhe  new 
techniques  which  have  been  developed  for  airplanes  to  the  response  of  air¬ 
ships  to  gusts.  The  gradual  penetration  of  thr  airship  into  the  gust  front 
will  he  considered  and  rhe  resulting  responses  and  loads  determined.  Recent¬ 
ly,  the  application  of  generalized  harmonic  analysis  to  the  gust  response  of 
airplanes  has  enjoyed  some  measure  of  success.  The  determination  of  the 
response  of  airs^ioi  ta_random  gust  dis’urbonces  appears  to  bejnore  oppro- 
pr:ate  since,  because  of  its  latge  relative  length,  the  airship  may  possibly 
encounter  a  succession  of  gusts  rather  than  a  single  "critical  gust". 

The  loads  acting  on  a  typical  airship  experiencing  discrete 
and  random  gust  disturbances  w  .1  be  determined  end  compared  to  the  load' 
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obioined  by  the  semi-empirical  approach.  On  the  basis  of  these  cat¬ 
enations,  critical  jdesign  parameters  for  the  envelope  and  empennage 
of  non-rigid  airships  rurrentiy  in  operation  will  be  derived. 
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CHAPTER  It 


EQUATIONS  OF  MOTION 


The  present  problem  Is  to  determine  the  verticol  ond  pitching 
motions  ond  the  structural  loodj  of  an  airship  resulting  from  a  vertical  gust 
disturbance  Wq  .  This  gust  is  considered  uniform  in  the  direction  normal 
to  the  plane  of  motion  of  the  airship.  It  is  Or, her  assumed  that  during  gust 
encounter  the  controls  of  the  airship  remain  lucked  and  its  forward  velocity 


V  is  con5tani. 
o 


The  coordinate  system  of  the  airship  in  a  disturbed  ran  'ition 
is  shown  in  Figure  I  where  h  is  the  vertical  displacement  o*  the  center  ot 
gravity  and  A  the  pitch  cngle,  defineJ  posn-ve  u*  .  I*'  d(x,|) 

be  *he  distributed  loading  acting  cn  the  oirship,  then,  for  equilibrium  in 
the  vertical  d;rect:on 


-  r 


J  f(X,t)J!x  =  o 

s  L  .  / 


'troJv 


and  for  quilibrium  of  moments  about  the  center  of  grovity 


/  (Xc^  -  Y.)  •{*  -  ° 

J  A  (1  j/ a 


The  distributed  looding  is  derived  .Turn  q  consideration  of  the  disturbance, 
mofio.1,  and  inertia  forces  acting  onthe  airship  and  may  bo  written  as 
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. .  -  - 


(2.1) 


(2.2) 


f“‘ =  7 ~  '  '*7  -  *] 

(2.3) 

where  '/J'X.  ond  "*  ore  the  dtstfiburea  oeroaynam<e"Torc.» 

due  to  motion  and  disturbance  end  m(x)  ii  the  distributed  matt  of  the  oir- 
ship.  Substitution  of  Eq.  (2.3)  Into  fas.  (2.  1)  and  (2.2)  result*  in  the 
following  statements  of  equilibrium: 


/  =L„ 

**/ 


•{ if-  ) 


U'x  =AL  +  V* 


Jy.  / 


(2,*) 


(2.  5) 


where  Mq  is  the  moss  of  the  airship  and  I  the  pitching  moment  of  ineitia 
about  the  center  of  gravity. 


2.  1  Aesodynomic  Force?  Due  to  Motion 
2. 1 .  1  Lift  Due  to  Motion 

The  aerodynamic  forces  ond  moments  acting  on  the  airship 
performing  unsteady  motion*  will  be  ■: '..-rived  by  linearized  slender-body 
theory.  Slender-body  theory  wa»  first  applied  by  Munk  (Ref.  t)  to  the 
calculation  of  the  force:  on  airship  hulls  or.d  later  extended  by  R.  T.  Jones 
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to  slender  wings  in  steady  flow,  Ref.  18.  hn  excellent  review  of  its  de..'|-,y>- 
ment  and  Applications  is  given  in  Refs.  19  and  20.  The  major  assumption  mode 
in  slender-body  theory  is  that  the  momentum  of  the  flow  in  a  plane  normal  to 
the  free  stream  is  the  same  c.s  if  this  flow  w*re  two-dimensional.  This  assump¬ 
tion  implies  that  the  variations  of  the  g*»me‘ricol  properties  of  the  body  in 
the  stream  direction  are  small  and,  also,  that  the  angle  of  attack  and  resultant 
motions  be  restricted  to  small  amplitudes..  The  force  per  unit  uxlol  diston-  e 
acting  cm  i  m«s  ,j  og..-.!  .  i!..i  <  S-* — *im*  'M*  nf  <-fwinqe  of  the  mom¬ 

entum  of  the  cross-flow.  Referring  to  Figure  1  ond  assuming  the  angle  of  attack, 
a(x,t)  ,  to  be  small,  the  velocity  of  the  cross-flow  at  section  A  -  A  1$ 
VQa(x,;)  ond  the  mc.mentum  is  yO  f  j  (x)  VQa(x  ,  t)  ,  whore  y1 
is  the  mass  density  of  the  *luld,  f  j  (x)  the  apparent  area  of  the  cross-sec'-on, 
and  the  free  stream  velocity.  Therefore,  the  force  per  unit  distance  act¬ 
ing  on  the  airship  because  of  Its  motion,  positive  upward,  it 


^ £*ri  _  [) 


[ f>  fj*) 


Within  the  limitations  of  slendor-body  theory,  the  substantial  derivative  iwi 

the  form  )  and  Eg.  (2.6)  becomes 

1  2  £  ,>i  /  —  .= 


7T  -(§*  «*•*!] 


where  the  angle  of  attack  at  each  cross-section  is 


ou<,i)  - 


0  -JL  d-  ( x-  yco  )  0 
V,  V. 
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Subttltutlng  Eq.  ^2.  8)  Into  Eq.(2.  7)  ond  dlfferentlotir.p,  the  lift  due  to 
motion  per  unit  length  b 


—  —  2 £  if  j  yt.  &  (*}  /-  (*-Yc0)  •ft (/>J 


...  A  J'(t)  _  £__  C  Ye  6-  -Y)  f,  (%} 
\/x  V02 


where  f  f*  Vi  ,  integrating  Eq.  (2.9)  clonq  the  length  of  the 

airship,  the  total  lift  due  to  motion  b 


=  2fi(6-AjV,U  +  £.[ t)T  ( 'yr  - 


/.<*)</*  -  A-  )/,(t)e/y-JL  ((*--*)£(«)*/* 

K/  7  W  " 


k>j 

/n 


(2. 10) 


In  applying  the  limit*  of  Integration  to  the  flat  and  second  term*  of  Eq.  (2.9) 
the  upper  limit  It  token  at  the  location  of  the  tectlon  of  moximum  tpan  of  the 
•oil  fin,  h  *  Xj  ,  because  the  turbulent  nature  of  the  flow  oft  of  tbit 
tectlon  dlvuph  the  twc-dlmemionai  Idealization  of  the  cmti-flow  (Ref.  18). 
The  quantity  f*  J fj(x}dx  it  recognized  ot  the  apparent  mots  of  the 
olahlp  and  mcy  be  written  ot  a  multiple  of  the  physical  mow,  Mq  , 
where  It  the  coefficient  of  additional  apparent  man  In  the  trontvers* 
erection.  The  quantity  -f(  ^  _  y)  £<*J  «/*  b  approximately  zero 
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bee  iure  f  j(x)  is  symmetrically  distributer  'bou*  the  center  of  gravity; 
hence  the  final  term  ir  Eq.  (2.  10)  will  be  omitted.  With  these  modificat  l-r 
Eq.  {2.  10)  may  be  rewritten  as 


-  2j_f{ (e~j_y/.)T  +6  (£)TfXr -*<*)] 


J 


v-  4  At  V,  e 


(2.11) 


where  (  f  ^ ^  is  the  value  of  the  apparent  area  of  the  cross-section  at 


T  • 


2. 1. 2  Pitching  Moment  Due  to  Motion 

Tie  pitching  moment  due~to  mot  I  an  about  the  center  of  grc\  itv 


A4  =  !  (ycq-  -x) 

Jkcjy  St 


(2.12) 


Substituting  Eq.  (2.9)  in  Eq.  (2.  12)  and  inieqroting  up  to  k  =  as  ir 
Eq.  (2.10) 
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l  =» 


/./**  T 


}  (Vo,) 


(?.  15) 


Z’  _  /Vg>fg/vr 


(?.  16) 


the  following  acfoc/nomic  ‘ability  derivative!  are  obtained-from  Eqi.  (2.  1 1) 
ond  (2.  .4): 


PC  = 

*  (l)-r 

77 Ip75 

j  V 

11 

Z  (.-r-Xtz!  V,)  r 

VoL 

PCf 

Z 

vOL. 

Cti  _ 

~  2  (X-r-  V'Cr)*  &)r 

c?Zr^J 

(\/t,) 
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^  „  .r  .  ^  _  rf.  .  , 

•«*,.  V*  ■\  ""v  •’/•"/.“.V. 


v/v,  v^v* %  v-  ^  , 


where 


(2.  ia 

Insteod  of  oiing  stobility  derJvoMves  os  expressed  by  &j».  (?.  17) 
above,  value*  based  on  current  design  practice  which  include  viscosity  effect* 
may  be  used.  Therefore,  o  more  rat  tonal  approach  ir.  expressing  the  lift  and 
moment  due  to  motion  at  small  angles  of  attack,  is  tc  wr'te  Eqs.  (2. 11)  and 
(2. 14)  os 


l/i  , 

(o—  o  (  Vol)  & 

v,/  r  jn»)  v, 

*  • 

+  Jz  \  VL  G  -Jx  J, 

Y/i  . 

/%  =,  C  (Vol)  SC*  (G  ~J.)  +  Q  (Vol)  <*C*  & 

*  -  vi /  r-  c/troi  v* 

-43Le" 

(2.1^ 

A  comparison  is  given  in  Table  I  of  the  stability  derivatives 
obtcir.od  from  E"t.  (2. 17)  and  by  semi-empirical  means  for  *he  example  air¬ 
ship  of  Chapter  IV.  The  theoretical  values  s^ceed  the  empirical. values  by 


Tfc 

•n 


«r  i 
\ 


‘  m}  w 


'<•**» 

‘V. 

/* 

s 

JM 

y 


approx imat ply  five  to  fifteen  percent,  which  may  be  attributed  to  the  neglect 
of  viscosity  in  slender  body  theory,  The  values  civen  '.<■  column  P  will  be 
used  in  this  report. 

2.  3  Aero  'ramie  Forces  Due  to  The  Gust  Disturbance 

The  lift  due  to  a  vertical  gust  disturbance  may  bt=  fuund  by 
consideiing  Wg(*,t)  os  the  cross-flow  velocity.  Replacing  V  Qa(  x  ,  t) 
by  W  q  ( x  ,  t )  in  Eq.  (2. 6) ,  the  lift  per  unit  distance  due  to  the  gust  dis¬ 
turbance  positive  upward,  is 


(?.  20) 

It  is  assumed  that  the  gust  profi'e  does  not  change  v.-ith  rime  but  travels  along 
the  x-direction  with  a  uniform  velocity  Vq  and  :$  derin«d  by  the  relai ion- 
ship 


%  <*>*>  - 

(-’■21) 

This  relationship  is  a  solution  of  the  general  wove  equation  and  has  the  prop¬ 
erty 


JL  %(V0t~*)~(±  +  V.  0_)  Wi(vA-y)  =  O 
O  t  Ox' 

r:  rv 


v ■, 

k 

*  i 
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which  simplifies  Eq.  (2.20)  to 


~  P  %.V0  -  Z  q  Wo  J1/.?*} 

**Jy  *  x  V.  w9  Jx 

where  Wq  is  the  amplitude  of  the  gust  profile.  Integrating  Eq.  (2.23),  the 
total  lift  Hue  to  the  gust  disturbance  Is 

L6  Z  9  f  We JjAi-xJ  J£(t)  Jy 

V‘  L,y  W‘  <** 

and  the  pitching  moment  about  the  center  of  grav.!/,  positive  nose  up,  is 


(2,23) 


(3.24) 


-  Zf  j  -  *)  (Y,t  -  x)  <//  lx)  t/x 

^  - 

In  terms  of  the  expressions  for  aricj  <VTC 

given  by  Eq.  (2  17),  Bqs.  (2. 24)  and  (2.  25)-moy  be  rewritten  as 


(2  25) 


-  v(v*l)  £j± 

r  A  14 


16  ,  <?  (Vol)  % 

"  '  Vo  -* 


(2.  26) 
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where 


£ 


%A M-*)  Ji,  Cxi 


#%. 


i aU 


Jy 


% 


/ 

sir 


i>oJ\ 


(*CZ  -  X)  ^  c/x 

U4  </x 


V  ’  (2.27) 

As  before,  more  realistic  values  for  and 


may  be  used  in  t:qs.  (2. 26)  In  place  of  the  values  obtained  by  slender-body 
theory,  Eqs.  (2. 17). 


2. 4  The  General  equations  of  Motion 

The  rquations  of  motion  of  the  cirship  in  distarbed  fl:gnt, 
fqs.  (2.  4)  and  (2. 5),  becone 


f  Vo  J(rp)  U  J 


-4x  Mo.  V*  0  =  q(Vol)  ^  Cl-  w,  </l 
f  J*  V.  * 


Tc  o<-4,)6  -9(v°l}[^1  (s-dL)  /  (Y.i)nd4L  / = 

r  •**  Jl*e)J 


(Vol)  Sfy  Wo 


$  l  VLJ  "  *T  vyo  Y>» 

*  s*  V,  * 


1 2, 28) 
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Th«*te  equations  may  be  simplified  by  introducii^  rb«  nondirrensiono!  time 
voriuble  si  -  which  ii  »he  oirshp  length  trnvellrd^and  the 

nondi.nenjional  vertical  displacement  J  r-/^.  .  In  transforming  to 
the  new  time  variable  si  ,  the  following  relationships  also  hold 


1 PA  V*  f  ' 

<4*  4-t  a* 


J 


e-  *r#' 


(?  29) 

wh«*re»{  )'  denotM  differentiation  with  respect  to  s)  Furthermore,  the 
phytlcal  and  aerodynamic  characteristic*  of  fbo  air, hip  may  bo  represented  by 
the  following  dimension!***  parumetets 


J-  D  ^<V*L)*/r 

2 '  **•  K  f'-*x) 

z  *  nZ*  f" 

z'  ^  z.<w,; 

S C  /i/ct  ff) 

, 

«rZL,  //rreJ 

SC  /J* 

</( £e» 

r(^7.) 

srJ  tr 

/ 

P.30) 
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.  V*  rfL  to" 


,'vMv'V 

to".  >  .* 


>  *Jr 


-z.vv 
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•n  feimj  of  tb«*  new  time  voriobla  si  ancLthe  (Kvomotefj  given  in 
Fq.  (2.  30)  obove,  ihe  equation!  of  motion  of  the  oi'Oiip  in  nondimenjional 
f  nrm  become 


The  qrowdr  of  lift  Hi.  ‘  to  a  sharp-ediei  gust,  GO  > 

ii  obtained  f-nm  Fq.  (2.3  )  by  putting  IV<» /-» -  f) -•  /  o  and  integrating 
Up  to  the  gust  front  /i  \ 

-  A 


kill 

,Jr)°  sr  (/J, 


oisi  i  fr 


A  J  f 


(2.  34) 


Thj?  the  !if{  uoring  on  on  air$nip  as  if  penetrate  o  sharp-edged  gutt  varies 
*ifh  the  apparent  area  cf  the  cross-section  at  the  gust  front.  Similarly,  the 
growth  of  moment  due  to  a  sharp-edged  gust,  ,  i*  obtained  from 

Eq.  (2.33)  . 


=  M.  [< -  ?)  •££*! 


(^x)-r 


off 


which,  after  integrating  by  parts,  yields 


(2.  35) 


_  MJ  or  /.  O 


’<r° 


o' J  '  F 


A  a  Fr 


(2.  d/ ' 

The  growths  of  lift  and  moment  due  to  cn  arbitrary  discrete  gust  profile, 
w*  '*>/w0  are  obtained  by  applying  Duhamei’s  integral  to  the  indicinl 
functions  of  Fqs..  (2.34)  and (2,  36): 
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w*r»)  # 

£  -  s-e< 


IV# 


(?.37) 


<2  w=  **/  %">+(  f  fyJd.1% 

-£T  *°  Jo  <fr£  w>  J  c 


IM • .  (4-r)cfr 

Cr  * 


(2.38) 


where  (/"'  is  o  dummy  variable  of  integration.  Knowing  the  forcing  function* 
Ad/  and  It  i*  then  possibla  to  jolve  the  coupled  equations 

of  motion,  Eqs.  (2,31),  for  the  responses  J(*)  ana  0(4)  .  Since  the 
airship  is  at  rejt  prior  to  disturbance,  th®  acrompanying  initial  conditions 
at  A  ~  O  are 


JJ (C)  *=£"(<>)-  0(°)  -  @  (o)  -  O 


(?.39) 


2.6  Random  Gust  Excitation 

The  response  of  the  airship  to  rardom  gust  excitations  may  be 
determined  by  applying  the  concepts  of  gen^.-Ured  harmonic  analysis  which 
are  discussed  in  Refs.  21  through  23.  With  the  application  of  these  concepts 
the  response  of  the  airship  can  be  described  in  term*  of  statistical  or  averago 
quantities,  the  most  Important  of  which  is  the  mean  square  value,  it  is 
assumed  that  the  gust  rurbulence  ’*  a  stationary  random  process  possessing 
the  proporties  of  homogeneity  ond  Isotropy,  'itafionarlfy  aod  homogeneity 
imply  that  the  statistical  characteristics  ',f  the  turbulence  aro  inva.  iont  with 
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time  ar.d  »{'•'>  -*,  re’.pectivjely.  Isotropy  implies, thoi  (he  statist i cal  proper  ti<*s 
of  the  turbulence  me  independent  of  the  airship  flight  poth. 

It  has  been  shown  (Ref.  21)  thot  the  mean  square  value  of  a  res¬ 
ponse  of  a  linear  system,  -  g-  *  ,  may  be  foun'1  from  the  following  relation¬ 

ship: 

tr> 

/%<"'>/*  (^y  /  tjs-  )  cfi  t*/” 

__  (?.  AO) 

where  H^(  iu  )  is  the  transfer  function  or  frequency  response  function  of 
the  system,  and  (  u)  is  the  power  specfraTdensity  of  the  disturbance, 
both  functions  of  the  harmonic  frequency  u.  The  output  C~ *  may  be  any 
response  of  the  system  such  as  velocity,  acceleration  or  streis.  The  transfer 
function  is  the  ratio  of  the  amplitude  of  the  system  response  to  the  amplitude 
of  a  sinusoidal  forcing  function. 

As  in  previous  theoretical  studies  of  o:r  turbulence  (Refs.  22 
and  23)  the  power  spectra!  density  of  the  vertical  oust  velocity  will  be  approxi¬ 
mated  by 


/■*"  /)Z 

§(~r)  =  iP  J  '  <  *  1  ^  / 


(2.  A  I) 


where  is  the  mean  sqi>a<-*  value  of  the  vertical  gust  velocity  end  ^ 

is  the  scale  of  turbulence,  which  can  be  considered  os  an  oporoxima‘e  measure 
uf  the  average  eddy  sire  in  the  turbulence.  The  transfer  function  of  the  air¬ 
ship  Is  found  by  subjecting  if  to  a  contimxrus  sinusoidal  gust  pattern  of  »he  form 


*4  fx,H=  c. 


i  ~r 


(?.42) 


Equations  (2.41)  and  (2.42)  may  be  rewritten  in  terms  of  the  reduced  frequency 

parameter  £  -  as 

t* 


C4,A)  = 


r) 


(2.43) 


where  .  The  steady  state  lift  end  moment  functions  ore  cb- 

tuined  by  substituting  Eq.  (2.44)  into  Eqs.  (2.32)  and  ( 2.33),  ~ 

V2r  («,*)= 


r.«) 


^  ^  </  =  /%  rJ)  c 


cJfd 


where 


O*) 


J  f*‘*F£*W 
//JW*  ✓jT 
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V  "J-7V  V7>"V7 


^  (■>-)=  r, 


(Ah. 


c~c*rCfrt-rl^hf) 


(2. 46b) 


The  displacements  T  and  9  wiltolso  be  periodic  and  may  be  written 


r  —  t '  urt  i  Asdl 

>  ~  >o  c  —  >o  e 


s  =  a.'*4 


where  To  ”■  ®0  are  the  amplitude  of  the  displacements  and  are 


complex  quantities.  Substituting  Eqs.  (2.  47)  and  (2. 45)  in  Eqs.  (2.31), 
the  equations  of  motion  In  matrix  notation  become 


~  » ...  tr^  n 

? 


-/-/  (/rx+tfj  .)*  r. 


ft*) 


- 

t  Vt 

-/-L  -t«\  a  a  £(4) 


The  transfer  functions  with  respect  to  the  vertical  and  pitch  accelerations 
ore  obtained  from  Eqs.  (2.48)  through  the  following  relationships 
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*V>  >  '*>'**■  ^  V***.."* 


(2.47) 


(2.48) 


A/  *  \Vt*  "w*  *•*  V*  ’%*'  %*  HJ"  *2*  **  **  *  *  -  *  «V*rM»re  XI 

•-  -  ^-v.v.v.v.v. 


■VAv.v, . 


;  .  .  «r.  it -v--v-'v ^-v’  r  v, 


The  bending  moment  may  be  expressed  in  nondimensionni  form 
by  the  following  coefficient  which  it  commonly  used  in  airship  design  (Ref.  R). 


0*1  — 


3*1 


£  Lt.(  \Zol  J 


*/> 


(3.4) 


Substituting  Eqs.  (2.9)  and  (2.  23)  into  Eq.  (3.  3)  and  converting  to  the  non- 
dimensional  variables  yf  and  ?  ,  the  bending  moment  coefficient 

at  station  due  to  a  gust  profile  j*V  (a)  becomes  ~ 


J*0*  -  J£_  W, 
(v.i)'n  v» 


A 

j Tr[' £r] 

b  ^ 

(tiL')l,  v  xz  -  £ 

1  ~./V.  1  1V»/V4_  W./n  7^ 


/t 

T 

^/vi 


/ 


*yvi_  2  “ 

^/vfc 

•  k4 


(3.5) 


where  h  the  indiclal  function  for  bending  moment 


m *//&- r)  <ff  ou if. 

^  J 
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Jo  tfif 


£  *  fo 
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A 

IJ=  / '  ff- r.)(f-f<>)  £(rW 


Cf 1 

T(  =  /  <r-  /v<r-  re,)^n)Jf 

Jo 


(3.7) 


The  shear  may  also  be  presented  in  coefficient  form  by  defining 
a  shear  coefficient  as  -.-  _ 


/  (V»V 


(3.8) 
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3.  ?  >he.ir  nod  .-iendlrvi  Moment  Due  to  Rnndcvm  ('-  .-  fv-itation 

The  opplii-otio"  cf  generalized  hain.-nic  orvolysis  enables  one 

fo  find  (>  f  m«in  square  v-lu**-,  of  the  sb»or  and  t . ^ng  moment  05  functions 

of  the  it'-'l*  of  NwboUn.  e.  The  procedure  entail-.  fast  the  determination  of 
♦he  tramfm  functions  wilh  inspect  fo  shear  and  h-ading  moment  ond  ‘hen  the 
Application  of  fo.  (2.40)  to  obtain  the  mean  squcre  v-olues. 


The  transfer  function  *ith  respect  to  shear  is  found  hy  specifying 
the  gust  pattern  oj 


w. 


fa-  ?) 


U  '7) 

in  Eg.  (2.  22)  and  the  displacements  £  and  8  ot 


It 

5o  <? 

u  - 

c44 

i 

in  Eqs.  (3.9).  The  resulting  steady-state  'hear  evolution  at  station 
becomes 


r 


(3.  13) 


%\r~ 


■)  cfiA 


(3.  1^) 


«  5 
fa’ 


"i 


■J,  ? 
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»;.*  7a>,>.v  '<• 


1‘  ••-■••'><•*,  rc.-rising  Fq.  (3.14) 


f  _  s-t 

O  -  >0  « 


('  15) 


where 


S',  -*f  f  [&  */£„] 


(3.  15) 


the  transfer  furmtlnn  with  respect  tc  the  shear  coefficient  it 


He  . 

s  f  (v*l)V>  V.  (V.l)1'} 


(3.  17) 


The  expressions  S  ^  and  S  are  the  reol  anrl  imoqinory  component  of 
the  termi  in  the  bracket*  of  Eq.  (3.  14)  and  are  piven  in  Appe'viix  A. 


Applying  Eq.  (2.40),  the  mean  square  value  of  the  shear  coeffi¬ 
cient  at  station  }'a  it 


4  <  /  [£,.;£] r±+*A'JAs4 

V°  Jo  (HtAf 


(3.  Ifl) 


In  a  similar  manner,  the  steady  state  amplitude  of  the  bending  lament  is 
found  to  be 
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\>r,  r,  c_ y  jr.ir,  •CKr.^.Jc,  <r ,  w  v  v  /.  V,’  fTPls.  f.  V  r.V.V'A  /.V 
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£-  FjY&'Iz)  -jj  ,:] 

if  j 


*  Yu  r  2:  +j  y73  *  t2 

’^/V.  L  ^ 


Rewriting  Eq.  (3. 19)  as 


hi  =  z?U*z-  [(*«)*  “  Mf 


the  transfer  function  with  respect  to  the  bending  rffoment  coefficient Ti 


and  the  mean  square  value  is 


(3.19) 


(3.70) 


(3.2') 


f**  =  771  /  [( **)*  +(**%.] 

”  '  yQ-t  \s*  /  *  - - - .  7 


V  Vs 


(3.  22) 


The  expressions  for  (8M)jj  and  (BM)  or®  fl’Vf>n  'n  Appendix  A. 
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3.  3  Correction*  to  Account  for  Viicoo*  E(fectj 

The  expression*  for  thecsheor  and  bending  moment  coefficVoM, 
Cqs.  (3. 9),  (3.  5),  (3.  ]  8)  ond  (3.  22),  oro  consistent  with  the  loeding  u*'  i 
to  derive  the  equation*  of  motion  (Eq.(2.  31)  J,  if  the  stability  derive1  ■  > 
ond  apparent  mas*  »erm*  of  Eq*.  (2.31)  are  derived  from  tlendsr-body  tlif-vy 
and  do  not  include  the  modification*  sug^ejted  in  Section  2.2  to  account 
fc  viicou*  effect*.  If  viscout  corrections  <r  d  in  the  equotiom  r f 
motion  they  should  also  be  included  in  the  $h„  ■>  »,  oendir  -  n  nt  *t  ex¬ 

pression*.  A  mean*  of  doing  this  It  to  multiply  the  appropriate  term*  of 
the  shear  and  bending  moment  equation*  by  the  mtio  of  the  empirical  value 
to  »“e  *lender-bod^  voluo.  The*e  ,“orre<-*io.w  •  i  «•*<■!'  t*‘-  1  zti  :'~~t 

they  may  be  applied  directly  to  the  Integral*  —  I . .  ,  the  indicia; 
function*  e  ,  ard  the  function*  aQ  ,  of 

Ea*.(A.  14)  and  y  ,  8^  ,  of  Eq*.  (A.  !A).  Table  II  preient*  a  tabu¬ 
lation  of  *ugge**ed  correction*  to  the  corrmp^. Jing  slender  body  value* 
which  hove  been  applied  to  the  example  •dtjhip  Qf  the  next  cheater  an'1 
may  also  be  used  fo:  other  airship*  of  jim’lor  design. 
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CHAPTER  IV 


DYNAM'C  RESPONSF  ANDSTRUCrURAL 
LOADS  OF  A  TYPICAL  AIRSHIP 


The  dynamic  response  and  resulting  structural  loadj  are  computed 
for  a  typical  non-rigid  airship  under  the  influence  of  discrete  and  random  gust 
excitations  at  an  altitude  of  5,000  feet.  The  parameter*  required  for  the  equ¬ 
ation*  of  motion  are  "  —  — 


?  =.6434 

p  =4.114 
r,  “  1.776 


r  2 


n 


=  -  3.318 
=  .7362 

=  . 2923 


and  were  evaluated  from  the  data  presented  In  lobi*  I!!.  It  I*  »een  that  the 
variation*  In  the*e  oaro meter*  are  very  wrwll  foMhe  four  major  type*  of  non- 
rigid  airship*  currently  Ir,  operation.  (Table  IV).  Furthermore,  *lnce  the 
operational  altitude  of  the  airship  It  probably  between  *ea  level  and  10,000 
teet,  variation*  tn  the  air  derolty  will  also  be  *mall  Thu*  the  parameter*  In 
Eq.  (4. 1)  above  and  an  altitude  of  5, 000  feet  will  be  considered  a*  typical 
of  all  current  non-rigid  airthlpt.  Therefore,  a  parametric  study  will  be  con¬ 
fined  to  varying  the  gust  dhtu-bance  parameter*.  The  data  for  Tobies  ill 
and  IV  were  taken  from  Reft.  16  and  27. 


4, 1  Response  to  Discrete  Gust  Excitation 

The  transient  response  of  the  example  airship  I*  obtained  by 
solving  Eq*.  (2.31)  with  the  prescribed  initial  condition*  of  Eq*.  (2.39). 


•4 


v  V  *V  *„ 


Vv-'v. 


>.V.  u  V  *v 


A  ;• b  v  ' <*X  v  v-viv1  <v.  l*wVvr  *r> st^t* 0*V  ^vA?^ oV^  A^«  <^vv^JvV^ 

j 

*-\ 


The  oroivsis  enfoi1!  first  the  determination  of  tl'~  •larp-edged  gust  fun  'ions, 
Eqs  (?.  34)  ond  (2.36),  the  specification*  of  o  pint  profile  *’V|Cr 
and  the  evaluation  of  ,v  H)  and  W*  6i,f  from  Eqs.  ("V  37) 

orvi  (2.3P). 

The  arowthj  of  lift  and  moment  due  to  a  ihorp-edqed  qioi  de¬ 
pend  nn  the  apparent  area  of  the  cross-section  of  the  airship  at  the  t,m!  front. 
For  wUiont  on  the  envelope  between  the  bow  or<d  fir*,  0  i  f  i  fr 
the  oprjrent  areo  i»  the  orea  of  the  cross-section 

fifr)  r=  IT  ■  o  £  f  £  fr 


(<  3 

where  R  i*  the  radius  of  the  envelope.  For  the  »~il  fin  »  onfiguratbn  of 
Figure  4a  the  apparent  orea  of  each  cross-section  i$  (?.cf.  24) 


f.i 

At 

.»■ 


The  fins  of  the  example  airship  also  have  a  crucS^ym  configuration  but  ore 
at  45  degrees  with  the  iortical  plane,  as  shown  in  Figure  3b.  SJnce  the  fins 
are  rotated  at  45  degrees,  the  vertical  lift  of  the  cross-section  shown  in  Fig- 
we  To  <i  Identical  to  t ho*  shown  in  Figure  3B”  if  bo'h  ore  immersed  in  the 
same  flaw,.  Hence,  Eq._(4.3)  will  be  used  in  tf>  >■  report.  The_ex  press  ion 
for  the  wTueopo  radius  P.  was  obtained  from  Ret.  25 
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Us  Ini;  Eq.  (4,4),  the  indicia!  function  *<x  lift  ond  momr-'* 
dye  to  tj"  penetration,  Fqs.  (2.34)  ond(2.36),  :i"' 


2 

(rl.t  -  \/.V/P7A  ~F.93/7^i'  +/  /f)tA3 

Mr  "  L 


./<?<?£  A 


2  r 

vrr^^l .  c 

(Al7  t 


</£  VTT/fV^r.  iVCiOA  ~/.V7.  of.  +2./ (.(>7 A 


-.  9'7$0  4  V.  /9^tA  rJ 


where 


ft  *■- 


J7F7  x/®  ft' 


For  cross-section*  on  the  fin  body  combination,  —  f  -  . 

Eq,  (4.  3)  wos  evaluated  from  F'guru  4  which  presents  the  approximate  plan- 
form  of  the  example  airship  fin.  The  apparent  area  in  this  region  was  approx¬ 
imated  by  straight  line  segments  as  shown  in  Figure  i  and  the  indicia!  funcfians 
for  this  region  ore 


LP  ,,,  .  /4oo  F  VS;/}/  (A-fr) 
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frF/  <0r 


(Mr 
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(4.9) 


where 


r_ 


*  .  Pt71 

J *T  "*  .  e*M 

Jr  *.93H 
^r*  ».4r«72 


(4.10) 


Plo»*  of  ond  ore  pr#>t«*nted  In  Figure  6. 
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The  discrete  gu;t  profile  selected  foi  this  study  is 


M-  f'i)  =  4-  r '-  c°-  * -~  i  o/r<  .^4 

tV,  Z  L  -  sig  ]  - 


Wi  MJ  =  0 


sS  ^  J  4 


(4.n) 


where  is  the  gust  gradient  distance  in  terms  of  the  airship  length. 

Plots  of  the  growths  of  lift  and  moment  due  to  this  profile  wero  obtained  by 
numerical  integration  of  Fqs.  <2. 37/  and  (2.33)  ond  arc  presented  In  Fig¬ 
ures  7  and  8  for  values  of  4*  equal  to  1/12,  1/A,  1/3,  1/2,  3/4. 


The  equations  of  motion  were  solved  numerically  by  using  cn 
Adam’s  method  expansion  retaining  third  differences.  The  required  first  four 
points  of  the  expansion  were  obtained  by  a  Runge-K"tta  solution  of  the  two 
differential  equations  of  motion.  This  numerical  procedure  Is  commonly  used 
in  solving  equations  of  this  nature  and  may  be  found  in  a  standard  text  such 
as  Ref.  26.  The  computations  were  performed  on  t^e  Burroughs  E-TOT  Digital 
Computer. 


A  typicai  time,  history  of  the  nondimensionol  response  of  the  air¬ 
ship  to  a  one-minus-cosine  gust  with  a  gradient  distance  -  1/2  is  pre¬ 
sented  in  Figures  9  and  10.  This  gradient  distance  was  selected  as  an  illus¬ 
tration  becom-*  it  fU.e  critical  loading  v/n  ihe  airship.  This  is 

discussed  in  the  next  chapter.  The  results  of  the  transient  response  ore 
summarized  in  Figure  11,  which  presents  the  maximum  values  of  the  non- 
dimensional  vertical  ond  pitching  accelerations  of  the  cenr<-  of  gravity  cs 
functions  of  the  perameter  .  The  times  for  these  accelerations 

to  '?sch  their  peak  values  are  also  presented  and  given  in  te/ms  of  the 
penetration  distance  y4 
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Ti*ri»  hirtorie«  cf  the  follow  -m  lwvp  olio  bw«n  rc.lrulot-' 
a  function  of  'bn  parameter  A j  : 


a. 

Figure  12. 

The  anal 

e  of  attack  ot  the  center  of  p- 

of  the  fin  ( 

Fcr 

-  .  8543)  colculoted  front 
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b. 

Figure  )3. 

The  non- 

dimensional  vertical  accelera? U* 
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i  cr 


(4.12) 


of  the  toil  jecTtori  al  the  center  oF'prejiure  of  the  fin, 
calculated  from 
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4.2  Response  to  Random  G ust  Excitation 

The  procedure  In  the  random  gust  analysis  requires  first  the 
determination  of  the  steady-state  lift  and  moment  functions,  F  j  (k)  and 
Fj(k)  »  9'v«n  by  Eqs.  (2.46)  and  the  computation  of  the  transfer  functions 
44r  and  /rj>  '  ,  Eqs.  (2.49),  Plots  of  thn  absolute  values  of 

/ F<‘*'/  ,  /K(A)/  ,  AV  7  1  ,  and  /W*  '/Z  are  given  in 

Figures  14  through  17  crs  functions  of  the  reduced  frequency  It.  The  mean 
square  values  J~"  and  6  "  ore  presented  In  Figure  18  as  functions  of 
the  scale  of  turbulence  ,  obtained  by  numerical  integration  of  Eq.  (2.50). 
The  integrals  v,  -e  found  to  converge  at  It  =  100.  The  power  spectral  density 
of  the  disturbance  used  for  the  computations,  Eq.  (2.43),  is  presented  in  Fig¬ 
ures  19  for  values  of  the  scale  of  turbulence  ranging  from  100  to  2000  feet. 

4.3  Structural  Loads  Due  to  Discrete  Gust  Excitation 

The  physical  and  cppnrent  mass  data  required  for  the  structural 
analysis  is  presented  in  Table  V  for  various  stations  J,  along  the  longi¬ 
tudinal  axis  of  the  airship.  The  dirtrlbuted  mass  of  the  example  airship  was 
taken  from  Ref.  16  and  is  presented  In  Figure  20. 

A  time  history  of  the  shear  coefficient,  Eq.  (3.9),  at  station 
T  o  =  .8171,  is  presented  in  Figure  21  as  a  function  of  si#  .  At  this 
station  the  shear  represents  the  normal  load  on  the  fall  section  and  is  of  primary 
importance  in  design  work.  A  time  history  of  th  »  bending  moment  coefficient, 
Eq.  (3.  5),  was  computed  at  y*  *»  .  4572  where  if  was  anticipated  to  be  a 
maximum  and  is  presented  in  Figure  22,  The  maximum  vo I ue  of  is 

very  important  to  the  designer  since  it  dictates  to  a  large  extent  the  design 
of  the  airship  envelope. 

The  following  computations  are  also  presented: 

a.  Figure  23,  Peak  valuator  /*'vV#  _at 

fo  *.4572  and  Cf /w./v»  at  =.8171 

os  a  function  of  *4$  ♦ 
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b.  Figures  24  and  25.  Time  his  ter  !oi  z,' 

5°  =.3171  ond  C»*/*»/V*  at  f0  •=  . 45/.' 

for  .d*  *  1/3,  to  Illust.-o*e  f!.«  relative  contributions 
-  of  the  dUtirbance,.. motion,  and  inertia  forces  to  the 
shear  and  bending  moment. 

c.  Figure  26.  First  positive  peak  ond  first  negc’<  <i  peak 

_  of  of  stations  f o  =.-2065,  .3540, 

,5900,  .7030,  .8171. 

d.  figure  27.  Envelope  of  pnifc  voluet  of  ^-9^/w*fv.  or •■! 
times  of  occurrence. 


4,4  Struct  ^ r!  loads  Due  tc  P.crdcrr.  GuitEU  citation 

The  numerical  results  of  the  rnn  ksm  gust  analjrsis  ore  summarized 
in  th a  following  figures: 

a.  Figure  28.  Square  of  absolute  value  6f  transfer  function 
with  respect  to  shear  coefficient  versus  Is  at  station 

=  .8171,  calculated  fra.t  feq.  (3.17). 

b.  Figure  29.  5quare-of  absolute  value  of  t-onsfnr  fund  Ion 
with  respect  to  bending  moment  coeffidert  versus  k  ot 
station  jTo  *  .4572,  calculated  from  Eq.  (3.21). 

c.  Figure  30.  Mean  square  values  of  shear  and  bending 
moment  coefficients  versus  scale  of  turbulence,  ealeu1- 
OTvw  i>  y  no  terlral  integration  of  Eqs.  (3. 18)  and  ( 3.271. 
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CHAPTER  V 


DISCUSSION  AND  CONCLUSIONS 

S.  1  Discussion  of  Col  ,i;loted  Reu/ts 

In  the  previous  chopter  the  rejponses  and  s>ructu'al  toads  of  the  example 
ai.  jh'p  experiencing  a  one- minus- cosine  gust  and  random  gust  disturbances  Hove  been 
calculated.  The/  are  presented  tn  nondi  mem  tonal  form  independent  of  the  gust  nmplt" 
tude  W  and  the  I'o'-word  velocity  of  the  airship  V  .  As  indicated  in  Chooter  IV. 

O  '  O  *  f 

there  results  also  rep r«*s*nt  approximately  the  characteristics  of  other  non-rta'd  ai?sh!ps 
currently  in  ope«-ntion. 

The  significant  observations  that  may  be  made  from  Figure  11  ore  that  l>'» 

pitching  acceleration  reaches  a  critical  value  for  a  pus!  .-ifh  a  pTodient  distance  of 

approximately  1/4  whe^-as  the  vertical  acceleration  is  critical  for  a  gradient  distance 

of  approximately  V4.  In  both  cases  the  peak  accelerations  <v  cur  shortly  after  the  pi/.* 

peak.  The  results  ot  the  random  gust  analysis  (Figure  1"  Indicate  that  a  scale  of 

turbulence  of  apprnwlmotnlv  1-1/2  airship  lengths  is  critical  for  the  vertical  accelarc’ion. 

Though  there  :j  insufficient  information  In  Figure  18  to  "’-.termine  a  critical  scale  of 

turbulence  for  the  pitching  acceleration.  It  appears  to  bn  test  than  1/2  on  airship  ierc*h 

Correlation  between  the  discrete  and  random  gust  results  it  favorable.  Both  analytet 

show  that  the  critical  wave  length  of  atmospheric  disturbance  is  shorter  for  the  pitching 

acceleration  then  for  the  verttcal  acceleration.  Alio.  the  values  of  ~  I  y^M/l  n,*> 
if"']  v 

j  agree  very  well,  the  ratio  of  the  former  lo  the  luftm-  being  1.5  to  1.3.  It 

should  be  noted  that  this  last  result  is  to  be  expected  a'»d  is  discussed  further  In  Re*.  ?2 

(r.  which  the  discrete  and  random  gust  alleviation  factors  of  a  rigid  airplane  were  f’Ua 

comparable. 

The  result*  of  the  structural  Jcxati  analysis  are  even  more  significant 
since  the  designer  is  primarily  interested  in  the  maximum  bending  moment  on 
the  envelope  and  the  maximum  load  on  the  toil  swcUoo.  The  results  in 


im- 


B' 


i  r 


Fibres  ??,  23  nod  27  shov.-  that  the  •rax  imum  b"i'd’ng  moment  on  rhe  eii^lon'* 
occurs  very  near  the  iamtion  of  the  center  pf  gravity  of  the  aifjhi.i  if*  MV, 


ar.-j  is  caused  Lv  a  gust  of  approximately  one  a'r«hip  length  (  A  &  -  ]/?). 

Figure  23  indicates  that  in  each  case  the  peak  of  the  bending  moment  o'  if  e 
center  of  gravity  occur*  shortly  after  tha  gust  peak.  The  norirol-  load  ..t  .'he 
tail  lection  (the  ihear  at  =  .5171)  olio  reochei  a  critical  value  for 

^  O  -  1/2  (Figures  2'  and  23).  This  critical  value  ii  attained  at  a  pene- 

trof  :on  distance  of  1. 26  airship  length*  or  ot  oppi  ox imately  half  an  oiuhip 
length  ottf-  *h«  oust  front  itrikei  the  tail  fin  leading  edge.  Aj  with  th» 
case  of  the  vortical  and  pitching  accelerations,  there  is  olio  favorable  ccrre- 
lct  ion  between  the  discrete  and  random  lead  analyse*.  The  random  analysis 
(Fiaure  30)  shows  that  es  scale  of  turbulence  of  approximately  one  amthip  length. 
(  /l  ~  1)  is  critical  for  both  the  maximum  erv/e'ope  bending  moment  ana 
the  toil  load.  Therefore,  both  meibodi  of  analysis  indicate  thct  the  critical 
wave  length  of  atmospheric  turbulence  Ii  the  some  for  the  envelope  and  tail 
The  mSAimum  values  of  the  curvet  or  Ficures  23  and  30  compare  os 


loadings. 
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The  results  of  tb;s  section  indicate  ttvjt  to  determine  a  critical 
gust  !'*natn  ;or  the  oirih'p  i'  is  necessary  to  go  beyond  the  determination  of 


the  response  of  th»  nirshin;  !.  e.  t  it  •;  necessary  fa  ca!-:;!at“  thn  loads  art  th® 
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airship.  Thi>  point  is  supported  by  the  fact  that  the  critical  values  of  the 
vertical  acceleration  of  the  center  of  gravity  and  the  pitching  acceleration 
are  not  caused  by  the  some  gust  length  which  results  in  critical  values  for  the 
loading  (  A  ~  1/2).  furthermore,  the  angle  of  attack  at  the  center  of 
pressure  r  f  the  fin  (figure  12)  indicates  that  A  a  =1/12  is  critical  whme- 
as  the  vertical  acceleration  at  this  station  (figure  13)  indicates  that  on  A a 
of  3/4  or  perhaps  greater  may  be  critical.  Therefore,  on  the  basis  of  the 
structural  load  analysis  it  is  concluded  that  the  critical  length  of  a  one-f'inus- 
cosine  oust  and  the  critical  scale  oi  turbulence  ore  both  approx irrotelv  one 
airship  length  and  that  Figure  23  could  form  the  basis  of  a  gust  design  criter¬ 
ion.  These  results  do  not  agree  with  Flomenhoft's  condition  (Ref.  15)  that 
a  gust  w'th  a  length  of  three  fin  chords  is  critical.  It  is  realized  that  his 
conclusions  are  based  on  the  acceleration  response  of  the  tail  instead  of  the 
loading. 


[Kijon  with  Sharp-Edged  Gust  Results 


A  comparison  will  be  mode  with  the  procedure  followed  in 
Refs.  16  and  17  in  v'hich-the  loadion  the  airship  are  analyzed  for  a  30  foot 
per  second  sharp-edged  gust  applied  norma!  to  the  piane  of  the  tail  fins. 

The  shear  and  bending  moment  distributions  are  co"-p.,fed  based  on  the  verti¬ 
cal  load  ot  the  stern  du';  to  this  sharp-edged  gust,  fhis  prccedurfltdoes  not 
consider  gust  loading  of  *he  envelope  and  presume-  that  this  stern  load  repre¬ 
sents  the  most  critical  gait  loading  scnditioriT  The  resuits  of  fh<T present 
Investigation,  however,  (Figure  22),  show  that  the  maximum  bending  mom¬ 
ent  occurs  before  the  one-mirus-cosin  e  gust  reaches  the  fin  leading  edge 
(at  A  =  .5//.  To  show  that  the  maximurvtail  load  as  computer1  by  the 
method  of  this  report  does  not  result  in  the  critical  bending  moment  on  the 
envelope  ond  that  the  envelope  loading  is  important,  the  fofowiag  calcul  ¬ 
ations  ore  made.  Assuming  that  WQ  ,  the  peak  of  the  one-mlnus-cosin* 
gust,  is  30  feet  per  second  and  VQ  ,  the  forward  velocity  is  122  feet  per 
second,  the  maximum  upward  tail  loed  is  found  from  Eq.  (5.  1)  to  be 


i 

The  moment  due  to  this  lood  about  the  center  of  gravity  is  1,440,00*0  ft -lbs. 
The  net  load  contributed  by  the  dlsirfbuted  loading  on  the  envelope  between 
the  center  of  gravity  and  the  tall  fins,  calculated  ot  the  Instant  that  the  tail 
lead  is  a  maximum,  was  found  io  bs  !4,62Q  lbs.  do~.-v.vsrd.  As  a  result  the 
net  bending  moment  ooout  the  centor  of  gravity  becom«  only  approximately 
200,000  ft-lb.  It  is  se^n  that  thfi  number  is  a  lrnol!  difference  between 
two  large  numbers  and  its  accuracy  depends  on  the  accuracy  of  the  distri¬ 
bute’  loading  on  the  rear  half  of  the  envelope.  This  distrlbui&d  loading  ts 
essentially  the  loading  due  to  the  gust  disturbance  which  itself  roquires  exper¬ 
imental  verification.  One  way  »o  accomplish  this  is  to  record  the  time  his¬ 
tory  of  the  bending  moment  at  the  center  of  prav!»v  and  to  verify  the  feet  *ha* 
the  gust  disturbance  at  the  tell  will  cause  very  high  local  shear  but  s  small 
bend’ng  moment  at  the  cvntw  of  gravity  station  of  the  alrshio. 

A  further  comparison  may  be  a f  the  maximum  bending 

moment  o«  the  envelope  and  the  shear  at  the  tali  section.  Reference  16 
gives 

(ff/*)***  =  g*°J  000 

5  =  -~/Oj  0OO  //s 

f 

and  the  results  of  the  present  investigation  give 

( 'e ■>/)»**  =  r-r-'ts 

//  C 
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fne  results  c(  S*  ?.,  16  then  ate  unconjervntive  in  comn^-ison  with  the  present 
investigation,  moreso  for  the  bending  moment  than  fcr  the  $hr>ar. 

In  calculating  the  operating  pressure  of  the  envelope,  the  value 
o»  (^et m  ) W  4 ft  that  is  commonly  used  is  .018  (Ref.  16)  and  is  inde*' 

.modent  of  the  gust  amplitude  and  the  airship  flight  velocity.  The  value  of 
((.0*}/> tjm  obtained  from  fcq.  (5.  1),  however,  is  a  function  of  both  the  gus1 
Linp! itude  and  the  flight  speed.  For  a  gust  amnlitude  o'  30  feet  per  second  and 
a  flight  soeed  of  172  feet  per  second,  Fp.  (5.  !)  gives  a  value  of  .0198. 


ft.  3  Can  etui  pas 

The  major  conclusions  of  this  report  me>  be  summarized  as 


foil 
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1. 


2. 


3. 


4. 
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There  is  no  apprecicb!e'd:ff;-ar'cc  bo^ween-lhe  discrete 
gust  factors  f  r  ,  ( '  .  and 

and  the  randvn  gust  factors”  ( 
and  ( (*»/"  /y.^ 

One-.ninirt-cosine  gusts  of  one  airship  length  end 
random  gust  disturbances  with  a  scale  of  turbulence  of 
one  airship  length  are  critical  for  both  the  envelope 
and  empennage  loading  of  the  uirship. 

Gust  loads  on  the  envelope  c-e  important  ond  should  be 
considered  in  a  theoretical  a^^lysis. 

The  critical  gust  lc*->d  Barometers  for  the  envelope  and 
empennage  of  rvon  rfgid  airship*  currently  in  operation, 
based  on  the  discrete  gust  analysis,  are 


= 


Ct  = 


43 


O  8/  W, 

XT 


30 


>4 


5.4  .■•»'dntior«  for  Futue  Research 

If  it  recommended  that  the  example  airship  be  instrumented 
for  flight  testing  to  obta.n  comoarison  with  the  (hoorerical  results  of  this 
report.  For  correlation  with  the  discrete  gulf  study,  time  histories  of  the 
followi'-g  should  bo  recorded: 

1.  #/«rfical  vslocify  and  acceleration  of  the  center  of 

gravity, 

2.  Pitching  angle,  velocity,  and  acceleration, 

3.  Vertical  acceleration  of  the  tail  section  , 

4.  Normal  load  on  the  empennage  and  bending  moment 
at  the  center  of  gravity  of  the  airship. 

For  a  correlation  with  the  theoretico1  rondon  gust  study,  the  transfer  functions 
of  the  vertical  ond  pitching  accelerations  arid  cf  the  empennage  load  and  en¬ 
velope  bending  moment  should  be  obtained.  During  rhe  re»t  Sur*  fh»  controls 
of  the  airship  should  be  locked  and  any  deviations  from  a  prescribed  setting 
should  be  recorded.  The  NACA  and  the  Boeina  Airpia.se  Company  hr.  « 
been  successful  in  conducting  similar  tests  for  airplanes  and  Refs.  28,  29  and 
30  are  recommended  for  further  study.  Two  basic  methods  which  hv.vt?  L.-,--n 
used  by  the  NACA  fRef.  28)  to  determine  the  o~ plitudes  of  the  transfer  func¬ 
tions  may  be  mentioned.  The  first  method  mdses  us*  of  il.a  power  ipecl.i.r...  of 
the  airplane  response  and  t!.e  gust  Input.  The  second  method  requires  the 
spectrum  of  airplane  response  and  the  cross  speefry-n  between  the  '  .cf  input 
and  the  airplane  response.  If  the  mousur'mnents  of  the  gust  inaut  ate  accurate, 
the  two  methods  should  give  Identical  results.  In  determining  a  means  by  which 
the  gust  input  could  be  measured,  consideration  should  be  given  to  the  fus¬ 
ibility  of  using  o  fiow~direction  vane  mounted  or>  o  boom  and  at  .’ached  to  the 
bow  of  the  airship.  Serious  boom  vibrational  problems  rr.oy  arise,  since  the 
boom  must  be  long  enough  so  that  the  vane  measurements  can  be  made  out  of 
the  airship  flow  field,  -’f  this  method  is  not-sucoesjful,  trie  method  suggested 
in  Rnf.  15  may  be  tried  in  which  an  inshumented  airplane  would  record  the 
turbulence  lave!  in  the  neighborhood  cf  the  airship  by  making  successive 
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poi\''<  '•  .■  fhc  oir;hip.  Tnis  lofter  method  i<;  feasible  Tot  measurement-,  of  th* 
po-.v"'  •-  ■  '•r.t'n  of  the  p-'*t  but  not  opplirable  far  ~'asun-tnei'ts  of  the  c 
spect'um  between  the  Oust  and  airship  ri’por,'.r.. 


It  is  realized  that  the  controls  locked  condition  is  somewhat 


idealistic  nd  that  the  pilot  is  constantly  using  the  controls  to  correct  for 
disturbance-..  I?  is  sygs»"*«fed,  therefore,  that  o  the-.-refica!  study  b*  mode 


which  will  include  art  automatic  control  system  ano  will  consider  the  effects 
of  control  surface  deflections  on  the  response  and  Icad.'g  of  the_airihlp. 


ASRl  TR  72-1 


* - 1  *V  -  -  -  -*  ’  *N  *•-  •  •-  ->  *> 


APPENDIX  A 


TRANSFER  FUNCTIONS  OF  THE  FX AMPLE  AIRSHIP 


A.  1  Verticcl  and  Pitching  Accelerations 

The  t.-apsfer  functions  with  respect  to  *he  vert:col  c.rxd  pitching 
acceleration*  of  the  an>Fip  ara  found  by  applying  Ctamer1*  rule  to  the  system 
of  Eqs.  (2.48): 


h  (A)  —  I  -  i  ( f>Y>  /  rvf  A  i  )J 

J. 1  . 

f:(A)  -/-  £  -  -c^axA 

_ F _ 

A  (A) 


*  } 

-  i A 


h  ( A  ) 


Fx  (A) 


Ark  J 

4-i  \  t\  f  ^ 


AU)r. 


-L\u 


A*  .  -> 

—  I  ~  - - A 

it  ^ 


The  lift  and  moment  functions  Fj(k)  and  are 
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Q]  --  1.4187 


a  j  ™  4.  8634 


o3-  3.5814 


a4  =  0.7232 


SubttituMna  Eq*.  (A.  4)  into  Eqj.  (2.46)  and  integiuim^, 
Fj(k)  and  F2(k)  can  b*  reduced  to  the  following  fofnn: 


/rcj)  _  /  y^jf/v  I  (A,  S/77  A  F.  !*  Coijff.  +4,  fW£ 

T  Mr 


A  /^v  S/s/  A  Ft  /  /&•  J  +■  l  f At  cos  A  />=*  -At  5//VA’F 


a 


di  cos  A  Fz  A  Av  soi  J  f-r  +  A  ■ 


/F  c A)  =  ^5  ^FFj-L  (&*  Fr~  A  CosAfp —  ^  CosAf* 

a  VU7 


A  A-  $///4f#  -&$  ccsAFr  a  Ft  >v.v4  -t*  8 ?) 

■f-  £  ^ fir  cos  A Fr  -  Ft  r/ssAF?  a  ^  s/*A F* 


i-  0+t  Cos  A  &f  A  &?S//V.Afr  +-8f,  cos  A  fr  C&g 
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T  \  '’ngularitle*  of  F  ^  ( k )  and  F k)  at  k  =  0  were  inv»tigat»d  by 
utHf/Sstg  the  jeriw  expomlan*  for  the  *lne  and  eoiine  term*  of  Eq*.  ( A.  *) 
and  (A.  7)  about  the  point  k  *  0  .  Fol losing  thi*  procedure 
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A.  2  Shoor  and  Binding  Moment 


The  real  and  imaciir>e»y  c^iponei.ti,  S  ^  and  S  .  m  ,  of  the 
steady-state  shear  coefficient,  Eq.  (3..  171,  ore  obtained  from  fcq.  f3.  14)  by 
introducing  the  displacement  amplitudes  o1  Eq.  (A.  1  1).  For  values  of  £ o 
between  0  and  fr  which  are  of  interest  in  this  irport: 
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Similarly,  for  the  jteady-statr-  Ijo-iding  moment  coefficient  , 
Eq.  (3.21),  $ub$titutiog  £q».  (A.  11)  into  Eq.  0- 19)  reiultt  in  the  following 
exprejjionj  for  the  interval  0  *  ~ 
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TABLE  I 

COMPARISON  OF  AERODYNAMIC  STABILITY 
DERIVATIVES  OF  EXAMPLE  A  IRSHIP 


A 
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Slender  Budy  Volue 
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TABLE  II 


VISCOUS  CORRFCTIONS  FOR  STMICTU--H  LOADS 


Shear 


Slender-B-xly  Term 
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C'<recMon  Slender  Tody  Term 
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AERODYNAMIC  AND  PHYSICAL  DATA 
OF  EXAMPLE  A 'POHIP 
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TABLE  !V 


AERODYNAMIC  AND  PHYSICAL  CHARAC?."  'STICS 
OF  FOUR  NON-RIGID  AftSHift 
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T6O0 

1590 

1620 

1750 

vr 
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-3. 32 

-2.  63 
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0.  307 
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0.1 71 
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0.00106 

0.00177 

0.00104 
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0.  732 
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1 . 30 
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TABS  E  V 

PHYSICAL  AND  APPARENT  MASS  •►-TEC-0  ALS  FOP 
SHEAR  AND  CCNDING  MOMFNi  EXPRESSIONS 
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FIG.  2  LOADING  DiA^AM  OF  AIRSHIP 
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FIG. 13  VERTICAL  ACCELERATION  AT  CENTER  QF  PRESSURE  OF 
FIN  D'JE  TO  I  “COS  GUST  _ 
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FI'?  14  ABSOLUTE  VALUE  OF  STEADY  STATE  LIFT  FUNCTION  F,  { 


FIG.  23  PEAK  VALUES  OF  SENDING  MOMENT  AND 
COEFFICIENTS 
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